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L-^ ■ The energy losses and spectra of Ultra High Energy Cosmic Rays (UHECR) 

are calculated for protons as primary particles. The attention is given to 
o ■ the energy losses due to electron-positron production in collisions with the 



microwave 2.73 K photons. The energy spectra are calculated for several 
models, which differ by production spectra and by source distribution, namely: 



^ ■ (i) Uniform distribution of the sources with steep generation spectra with 



indices 2.4 - 2.7, with cosmological evolution and without it. In this case 
it is possible to fit the shape of the observed spectrum up to 8 • 10^^ eV, 
but the required CR emissivity is too high and the GZK cutoff is present. 
(ii) Uniform distribution of the sources with flat generation spectrum dE/E"^ 
which is relevant to GRBs. The calculated spectrum is in disagreement with 
the observed one. The agreement at i? ^ 8 ■ 10^^ eV can be reached using a 
complex generation spectrum, but the required CR emissivity is three orders 
of magnitude higher than that of GRBs, and the predicted spectrum has the 



GZK cutoff, (in) The case of local enhancement within region of size 10 - 30 
Mpc with overdensity given by factor 3 - 30. The overdensity larger than 10 
is needed to eliminate the GZK cutoff. 

I. INTRODUCTION 

The energy losses of UHE protons in extragalactic space are caused by interaction with 
microwave radiation. The contribution of IR and optical radiation is small (for a detailed 
review of energy losses and the resulting spectrum see |]T|). The main contribution to energy 
losses is given by expansion of the universe, electron-positron pair production and pion pro- 
duction. The latter process results in steepening of the proton spectrum referred to as the 
Greisen-Zatsepin-Kuzmin (GZK) cutoff [§. The GZK cutoff is not seen in the observational 
data (for a recent review see [§]). The most conservative approach to explanation of observa- 
tions is astrophysical one: the protons are accelerated in extragalactic astrophysical sources 
(normal galaxies, compact objects in normal galaxies, e.g. GRB engines, AGN etc) and 
propagate towards us. This approach comprises three aspects: acceleration to UHE, total 
energy release in a source and propagation in extragalactic space. This most conservative 
approach is considered as (almost) excluded, with certain caveats, however. The models in 
which the GZK cutoff is absent or ameliorated include nearby one-source model (see [^Q and 
most recent work |^); the Local Supercluster model, in which the density of UHECR sources 
is locally enhanced ( [|T],01, for a recent work see |^); and finally widely discussed GRB model 
which, according to calculations 0, gives a reasonable agreement with observations. 

In this paper we shall analyse the two latter models. 



II. ENERGY LOSSES 

We present here the accurate calculations for pair production, p + 'Jcmbr — ^ p + e+ + e^, 
and for pion production p + 'Jcmbr ^ N + pions, where 'Jcmbr is a microwave photon. 

The energy losses of UHE proton per unit time due to its interaction with low energy 
photons is given by 

1 dE c /"^ , , . „, . f°° , nfe 



Edt 2T^J.J''''^''^f^''^'^U/'^^ (') 

where F is the Lorentz factor of the proton, e,- is the energy of background photon in the 
system where the proton is at rest, eth is the threshold of the considered reaction in the rest 
system of the proton, cr(er) is the cross-section, /(e^,) is the mean fraction of energy lost by 
the proton in one p7 collision in the laboratory system, e is the energy of the background 
photon in the lab system, and n{e) is the density of background photons. 
For the CMBR with temperature T Eq.([^) is simplified 
1 dE cT 
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From Eqs.(|l]) and (§) one can see that the mean fraction of energy lost by the proton in lab 
system in one collision, /(e^) = (1 — x) = {Ep — Ep)/Ep, is the basic quantity needed for 
calculations of energy losses. The threshold values of these quantities are well known: 

/pair ~ ) Jpion ~ i , n / ' \"^/ 

rup nip 1 + ler/mp 

where /pair and /pion are the threshold fractions for p + 7 -^> p + e+ + e^ and p + 7 -^ A^ + tt, 
respectively, and /x is the pion mass. 

For the accurate calculations of energy losses the fraction / properly averaged over dif- 
ferential cross-section is needed. 

Pair production loss has been previously discussed in many papers. All authors directly 
or indirectly have followed the standard approach of Ref. [|10| where the first Born approxi- 
mation of the Bethe-Heitler cross-section with proton mass mp —>■ 00 was used. In contrast to 



Ref. JTOl, we are using the first Born approximation approach of Ref. [|Tl|], which takes into ac- 



count the finite proton mass. We also use the exact non-relativistic differential cross-section 



from 1 12]. This allowed us to calculate the average fraction of energy lost by the proton in 
lab system by performing fourfold integration over invariant mass of electron-positron pair 
Mxi over an angle between incident and scattered proton, and polar and azimuthal angles 
of an electron in the cm system of the pair (see Appendix A for further details). 

Calculating photopion energy loss we followed the method of papers [l^Jl^ . Total cross- 
sections were taken according to Ref. ||15|. At low cm. energy E^. we considered the binary 
reactions p + 7 -^ n + N (including the resonance p + 7-^A), p + 'j ^ n^ + A^+, 
and p + 7 -^ p° + p. Differential cross-sections of binary processes at small energies were 
taken from |TH]. At Ec > 4.3 GeV we assumed the scaling behaviour of differential cross- 



sections. These were taken from Ref. |]T^. In the intermediate energy range we used an 
interpolation approach which allows us to describe the residual part of the total cross-section. 
The corresponding differential cross-sections coincide with low-energy binary description and 
high-energy scaling distribution and have a smooth transition between these two regimes in 
the intermediate region. 
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FIG. 1. a) UHE proton energy losses E~^dE/dt (present work: curve 1; Berezinsky and Grig- 



orieva (1988) f^^: curve 2; Stanev et al 2000 \lq] : black squares). The line 3 gives energy losses 
due to redshift (Hq = 65 km/sMpc). h) The derivative dhQ[E)/dE, where bQ{E) = dE/dt at present 
epoch z = 0. 

The results of our calculations are presented in Fig.^ in terms of relative energy losses per unit 
time E'^dE/dt as function of energy (curve 1). Also plotted is the derivative dbo{E) / dE , 
where bo{E) = dE/dt (Fig. |l|b). This quantity is needed for calculation of differential energy 
spectrum (see section III). In Fig. |l] we plot for comparison the energy losses as calculated 
by Berezinsky and Grigorieva 1988 [T^ (dashed curve 2). The difference in energy losses 
due to pion production is very small, not exceeding 5% in the energy region relevant for 
comparison with experimental data(-E < 10^^ eV). The difference with energy losses due 
to pair production is larger and reaches maximal value 15%. The results of calculations by 
Stanev et al |jl8| are shown by black squares. These authors have performed the detailed 
calculations for both aforementioned processes, though their approach is somewhat different 
from ours, especially for photopion process. Our energy losses are practically indistinguish- 



able from ||T^ for pair production and low energy pion production, and differ by 15-20% for 
pion production at highest energies (see Fig. |l]). Stanev et al claimed that energy losses due 
to pair production is underestimated by Berezinsky and Grigorieva |jT9| by 30-40%. Com- 
parison of data files of Stanev et al and Berezinsky and Grigorieva (see also Fig.|I|) shows 
that this difference is significantly less. Most probably, Stanev et al scanned inaccurately 



the data from the journal version of the paper |19 



III. UNIFORM DISTRIBUTION OF UHECR SOURCES AND GZK CUTOFF 

The GZK cutoff is a model-dependent feature of the spectrum, e.g. the GZK cutoff for a 
single source depends on the distance to the source. A common convention is that the GZK 
cutoff is defined for diffuse flux from the sources uniformly distributed over the universe. 
In this case one can give two definitions of the GZK cutoff. In the first one the cutoff is 
determined as the energy, Egzk ~ 3 x 10^^ eV, where the steep increase in the energy losses 
starts (see Fig. |l|). The GZK cutoff starts at this energy. The corresponding pathlength 
of a proton is Rgzk ~ {E^^dE/dt)'^ ?^ 1.3 • 10^ Mpc. The advantage of this definition of 
the cutoff energy is independence on spectrum index, but this energy is too low to judge 
about presence or absence of the cutoff in the measured spectrum. More practical definition 
is -£"1/2, where the flux with cutoff becomes lower by factor 2 than power-law extrapolation. 
This definition is convenient to use for the integral spectrum, which is better approximated 
by power-law function, than the differential one. In Fig.|^ the function E^^^^\j(> E) , where 
J(> E) is calculated integral diffuse spectrum, is plotted as function of energy. Note, that 
7 > 7g is an effective index of power-law approximation of the spectrum modified by energy 
losses. For wide range of generation indices 2.1 < 7^ < 2.7 the cutoff energy is the same, 
Ei/2 ~ 5.3 ■ 10^^ eV. The corresponding proton pathlength is R1/2 ~ 800 Mpc. 
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FIG. 2. Integral UHECR spectra with indices of generation spectra jg = 2.1 and jg = 2.7 (solid 
and dashed line, respectively). The vertical dotted line shows the energy E1/2, where the calculated 
flux becomes two times lower than power-law extrapolation. 

Using energy losses given in Section 2, we calculated the diffuse spectra for the model 
when sources are distributed uniformly in the universe. We followed the method of calcula- 



tion suggested in Ref. [|T9[. We use two assumptions for uniform distribution of the sources: 
(i) with evolution of the sources described by factor (1 + z)"^ in comoving frame [Q, and 
(a) without evolution. The power-law generation spectrum with index 7^ was assumed. 
We made different assumptions about maximum energy in the generation spectrum, namely 
-E'max = 3 ■ 10^° eV , -Emax = 1-10^^ cV, and -Emax = C)0. Varying parameters •jg and m we 
fit the AG ASA and Akeno data [g^. 

The fit of UHECR data with help of evolving sources was made in the past (e.g. see 
Ref. [^] and |JI|). The widely used fit for the AGASA data with jg = 2.3 and with assumed 
mixed composition of galactic and extragalactic UHECR was found by Yoshida and Teshima 



2^ . Recently Scully and Stecker ||23| made calculations similar to that above for UHECR 
produced by GRBs. 



We calculate spectra using the formalism of Ref. [IS 



where Zg is a redshift at generation and Eg{zg) is energy of a proton at generation, if at 
present {z = 0) its energy is E: Eg{zg) = X{E,Zg)E and X{E,Zg) is calculated numerically 
using energy losses dE/dt accounted for their time evolution; Cq = n^Lp is CR emissivity at 
2; = (no and Lp are space density of the sources and their CR luminosity, respectively). As 
the general case we assume cosmological evolution of the sources given by C{z) = Co{l + z)"^, 
where the absence of evolution corresponds torn = 0. All energies in Eq.(^) are given in GeV 
and luminosities in GeV/s. Dilation of energy interval is given by [|T^ (see also Appendix 
B): 



^ = (l + %)exp 



d^ 1/2 fdbo{E') 



^0 V dE' ) E'=(l+z)E,(z,)^ 



(5) 



where bQ{E) = dE/dt is energy loss due to interaction with CMBR photons at 2; = 
(adiabatic energy loss due to redshift must not be included!). Derivative dbo{E)/dE at 
z = is given in Fig.([l|b). 

For particles with energies E ^ 1 x 10^^ eV the maximum redshift for evolution of 
CR sources 2;max is not important if it is larger than 4. Integration over large z gives small 
contribution when the generation energy Eg{E, z) reaches the value Eeq{zm), at which energy 
losses due to pair production and redshift are equal. Then the maximum redshift Zm{E) of 
the epochs contributing to the flux of protons with energy E is determined by equation 
(1 + Zm)E = Ecciizm)- FoT energies E < 1 x 10^''' eV the maximum redshift of the source 
evolution Zmax might be important. In these cases we fix it as Zmax = 4. 



25 - 




1 I I I I i 



17 



18 



19 



20 



21 



FIG. 3. UHECR spectrum as observed in Akeno (triangles) and AG AS A (filled circles) experi- 
ments. The curves show the predicted differential spectra for the uniform distribution of sources with 
or without evolution. The case without evolution (m = 0, jg = 2.7) is given by curves (1),(2),(3) 
for maximum generation energy Emax = 3 • lO^'^ eV, 1 • 10^^ eV and oo, respectively. The dashed 
curve 4 describes the evolutionary model with m = 4, 7g = 2.45 and -Emax = oo. 

We can fit the Akeno-AGASA data in both cases, with and without evolution. The 
spectra without evolution, m = can fit the data starting from relatively high energy 
E > 1- 10^^ eV. The fit needs 7^ = 2.7. The curves 1, 2 and 3 in Fig.^ show the spectra with 
different Emax equal to 3- 10^° eV, 1 ■ 10^^ eV and 00, respectively. The fit without evolution 
(curves 1, 2, 3) needs Cq = 4.7 ■ 10^^ erg/Mpc^yr, while the fit for evolutionary case (curve 
4) needs Cq = 1.3 ■ 10^^ erg/Mpc'^yr. The difference between these two emissivities is caused 
mainly by flatter generation spectrum in the evolutionary case. 

The required emissivities can be compared with most powerful local emissivity given by 
Seyfert galaxies Csy = nsyLsy Using the space density of Seyfert galaxies risy ~ 10^^^ cm^^ 



9 



and the luminosity Lsy ~ 10^*^ erg/s one obtains Csy ~ 1 ■ 10'^^ erg/Mpc^yr, which is almost 
4 orders of magnitude less than CR emissivity needed in no-evolutionary case and one order 
of magnitude less than one in the evolutionary case. 

As Fig.|] shows the models with uniform distribution of the sources are excluded by 
absence of GZK cutoff in the observations. They give a good fit to the lower energy data. 
This fit needs large 7^, and thus too large energy output of the sources, riL. It is possible to 
overcome this difficulty using an assumption that production spectrum is flat at low energies 
and has a steepening at some high energy E^.. Assuming, for example, that spectrum is 
oc E^"^ at low energy, and oc E^"^'"^ a.t E > Ec = 1 x 10^ GeV, one obtains the required 
CR emissivity ( see Section V) C = 3.7 x 10^^ erg/Mpc^yr, i.e. less than observed total 
emissivity due to the Seyfert galaxies. A plausible assumption is that the population of 
UHECR sources is comprised by galaxies with moderate activity of AGN, which at higher 
luminosities are linked to Seyfert galaxies and BL Lacertae. There were recently found the 
observational indications that the latter galaxies could be the sources of observed UHECR 



2J]. If such sources had large local overdensity, the GZK cutoff would be less noticeable. 



We shall study this possibihty in the next Section. 

IV. LOCAL OVERDENSITY OF UHECR SOURCES 

Local overdensity of UHECR sources makes the GZK cutoff less sharp or eliminates it 
[T|]. Clustering of galaxies is a gravitational property, which is determined by mass and not 
by internal activity of an object. The galaxies of the same masses with active galactic nuclei 
or without them, with burst of star formation or in quiet phase, are clustering in the same 
way. Therefore the optical catalogues give a reasonable indication to expected clustering of 
UHECR sources. 
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FIG. 4. The effect of source overdensity on UHECR spectra for Roverd = 30 Mpc and different 
overdensity ratios n/n^ = 1, 2, 10, 30 (curves 1, 2, 3 and 4, respectively). 

The nearby structure that can affect the GZK cutoff is Local Supercluster (LS) of galaxies, 
which has a form of ellipsoid with semi-axes 20 and 30 Mpc. The LS overdensity of galaxies 



is estimated by factor ~ 2 ( see p5| and references therein). Such overdensity does not solve 
the problem of GZK cutoff [0,^1. We shall calculate here UHECR spectra for different local 
overdensities ti/uq, where no is mean extragalactic density of UHECR sources. We use the 
various sizes of overdensity region Roverd, equal to 10, 20 and 30 Mpc. The results of our 
calculations are presented in Fig.§ for 7^ = 2.7, -m = and for four values of overdensity 
n/rio equal to 1, 2, 10 and 30, assuming the size of overdensity region 30 Mpc (the results for 
Roverd = 20 Mpc are not much different). From Fig.|^ one can see that overdensity u/uq > 10 
is needed to reconcile well the calculations with observational data. 



11 



V. UHECR FROM GRB 



In GRBs the protons can be accelerated to Ultra High Energies [^,0. The strong 



indication that UHECR can be produced by GRBs, the authors of Ref. p6|j2^ j9|] see in the 



equal emissivity £ in GRBs and UHECRs. Scully and Stecker [^] argue that in fact the 
energy output in cosmic rays is higher than in GRBs. We shall analyse here the problem of 
energy output combined with the spectrum shape. 

For energetically most favourable CR generation spectrum dE/E"^, advocated in [ PB| , P7[ , 
the diffuse spectrum of UHECR can be found as 

JpiE) = ^ -^ ^E-' r^^ dzJl + z,r-'l'X-\E, z,)^. (6) 

P^ ' 47r inf^ Hq Jo 9K aJ ^ ' ^' dE ^ ' 

The calculated spectra for non-evolutionary case m = and for evolution of GRB sources 
with m = 4 are displayed in Fig.^ by curves 1 and 3, respectively. The required CR emissivity 
is Co = 2.0 ■ 10'^^ erg/Mpc^yr for both cases. It is two orders of magnitude larger than that 
observed in GRBs [Hf Sgrb = 1 • 10"^^ erg/Mpc^yr. 



Apart from the problem of too large energy output, these models do not fit the ob- 
served spectrum shape and predict the standard GZK cutoff. To obtain the agreement with 
spectrum shape one can use an artificial E""^ spectrum with steepening at energy E^.. 

At energy E > E^. the generation spectrum of a source is 

Q,{E,^z) = ^^;^^^E^^^^ (7) 

while at E < Ec this spectrum is assumed to have 1/E'^ shape. It is easy to verify that this 
spectrum is correctly normalized to the luminosity Lp. The diffuse spectrum can be readily 
calculated at E > E^ as 

M^) - -^ ^» ilT^T^ L "-'^ ^ s)'"-=''A-^.(E. ./-^. (8) 

The fluxes given by Eq.(|D are displayed in Fig.^by curves 2 and 4 for cases m = and m = 4, 
respectively. These spectra agree well with the Akeno-AGASA observations at energies lower 

12 



than the GZK cutoff, but require higher emissivity, Cq = 3.7 ■ 10^^ erg/Mpc^yr (m = case, 
curve 2) and Cq = 3.1 ■ 10^^ erg/Mpc^yr {m = 4 case, curve 4). 

We conclude thus that UHECR from GRBs exhibit the standard GZK cutoff and require 
CR emissivity 2-3 orders of magnitude higher than that observed in GRBs. Our conclusions 



agree with that of Ref. [23 




log.oE, eV 

FIG. 5. UHECR spectra from GRBs for E~^ generation spectrum (curve 1 for non- evolutionary 
m = case and curve 3 for evolutionary m = A case) and for E^"^ spectrum with steepening (curve 
2 form = 0, E^ = I ■ 10^ GeV, jg = 2.7, and curve 4 for m = A, E^ = 1 ■ W^ GeV, 7^ = 2.45; 

VI. CONCLUSIONS 

We have performed accurate calculations of energy losses of UHE protons due to electron- 
positron pair production and pion production in collisions with microwave photons. The 
diffuse spectra of UHE protons have been calculated for uniform distribution of the sources 
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in the universe for different maximum energies of the generation spectrum. The generation 
spectrum with index 7^ = 2.7 provides a good fit for energy range 1 ■ 10^^ — 8 ■ 10^^ eV in 
case of absence of evolution. For the case of evolution the good fit is given by m = 4 and 
7g = 2.45 in the energy range 1 ■ 10^"^ — 8 • 10^^ eV. 

Local overdensity of UHECR sources, e.g. in Local Supercluster, can reconcile the weak 
GZK cutoff with UHECR data only if overdensity is larger than 10. The existing astronomical 
data favour much smaller overdensity, of order of 2. 

The excellent fit of observational data in energy range 1 ■ 10^^ — 8 • 10^^ eV by 
astrophysical spectrum is rather impressive. It requires too high local CR emissivity 
£0 ~ 1 ■ 10^^ erg/Mpc^yr (in case of an evolutionary model), but this emissivity can be 
drastically reduced down to 3 ■ 10^^ erg/Mpc^yr by assumption of E'"^ spectrum with steep- 
ening at ii^c = 1 ■ 10® GeV. Normal galaxies with weak AGN {quasi- seyferts) meet this energy 
requirement. The local overdensity of these galaxies could ameliorate the GZK problem, but 
required overdensity, n/no > 10, is larger than that observed n/no ~ 2. 

UHECR from GRBs have a standard GZK cutoff. The predicted E'"^ generation spec- 
trum gives bad fit to the observed spectrum at energy lower than GZK cutoff. In case the 
generation spectrum is modified to give a reasonable fit, the required CR emissivity exceeds 
that observed in GRBs by three orders of magnitude. 
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APPENDIX A: CALCULATIONS OF ENERGY LOSSES 

Pair production energy loss of ultrahigh-energy protons in low-energy photon gas, e.g. 
CMBR, 

P + IcMBR ^ p + e+ + e" (Al) 

has been previously discussed in many papers. The differential cross-section for this process 



in the first Born approximation was originally calculated in 1934 by Bethe and Heitler [29 



and Racah [^. In 1948 Feenberg and Primakoff ||31| obtained the pair production energy 
loss rate using the extreme relativistic approximation for the differential cross-section. And 
in 1970 the accurate calculation was performed by Blumenthal [|10[. Later some analytical 
approximations to differential cross-sections were applied to this problem in Ref. ||32|| . 

All authors neglected the recoil energy of proton putting mp -^ oo, the effect being 
suppressed by a factor of me/nip ^ 5 x 10"'*. 

In spite of the fact that all calculations actually used the same Blumenthal approach, 
there are noticeable discrepancies in the results of different authors; they were clearly demon- 
strated in Fig. lb of the Ref. pf. 

To clarify the situation we recalculated the pair production energy loss of high-energy 
proton in the low-energy photon gas. In contrast to Ref. |T^ we use the first Born ap- 



proximation approach of Ref. |^ taking into account the finite proton mass. The exact 
non-relativistic threshold formula with corrections to different Coulomb interactions of elec- 



tron and positron with the proton (see e.g. Ref. [Q) was used. No series expansions of 
cr{E^) were involved in our calculations. 

Our strategy was to calculate the average energy transfer x = E' /Ep, where Ep and E' 
are the incident and final proton energies respectively, in the laboratory system by performing 
the direct fourfold integration of the exact matrix element over the phase space. It should 
be noted, that direct numerical integration, especially at high energies, is difficult in this 
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case because of forward-backward spikes in the electron-positron angular distributions. To 
overcome this problem, we performed two integrations over polar and azimuth angles in the 
e"'"e~ subsystem analytically. This was facilitated by using of the MATHEMATICA 4 code. 
The residual two integrations over energy and scattering angle in the initial pj subsystem 
were carried out numerically. We calculate simultaneously the total cross-section for pair 
production. The accuracy of our calculations was thus controlled by comparison of calculated 
total cross-section with the well-known Bethe-Heitler cross-section. 

The average fraction of proton energy lost in one collision with a photon is plotted in 
Fig. P as a function of the photon energy in the proton rest system. 
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FIG. 6. The average fraction of the incident proton energy Ep carried away by e~^e~ pair as a 
function of the photon energy E^ in the proton rest system. 



The product of this fraction and the total cross-section for pair production is shown in 
Fig. 1^. This function should be integrated over the photon spectrum to obtain the average 
energy loss due to pair production in the photon gas with this spectrum. 
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FIG. 7. The product of fraction of energy lost, (l—x), and the cross-section for pair production, 
'jp — > e~^e~p, or photopion production p + j ^ X as function of photon energy in the proton rest 
system E^. 



The comparison of our calculations with Ref. |]18| shows the negligible difference (see 

Fig.HD. 



APPENDIX B: CONNECTION BETWEEN ENERGY INTERVALS AT EPOCHS 

OF PRODUCTION AND OBSERVATION 



If we consider the protons with energy E in the interval dE at the epoch with redshift 
z = 0, what will be the corresponding interval of generation dEg at epoch z, when energy of 
a proton was Eg{z)l The connection between these two intervals is given by Eq.(36) of Ref. 
|[T9|1 . Here we shall confirm this formula using a different, more simple derivation. Note, that 



intermediate formulae (40) and (41) used for derivation of final Eq.(36) in Ref. |T^ have a 
misprint: the correct power of (1 + z) term there is 3, not 2. 

Regarding the energy losses of a proton on CMBR, dE/dt and {l/E)dE/dt, at arbitrary 
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epoch with redshift z^ we shall use the following notation 



W^^/cmbr' ' \E dt ^ 



(Bl) 



CMBR 



As it is readily seen from Eq.(0), the energy losses at arbitrary epoch z can be found as 



P{E,z) = {l + z)'Po{{l + z)E) 



(B2) 



where notation bolE) and (3o{E) are used here and henceforth for energy losses at 2; = 0. 
The energy loss due to redshift at the epoch z is given by 



if) „-^o(l--P- 

V / r— sh 



(B3) 



The energy of a particle at epoch z, 



dt 



'dE\ fdE" 

lit] , \~dt 



r— sh 



CMBR. 



E,{z)=E + j^ 
can be easily rearranged as 

with help of Eq.( |B2[ ) and expression dt = dz/ Hq{1 + z)^'"^ for time interval. 



(B4) 



(B5) 



Differentiating Eg. ( |B5D over E, one finds for energy interval dilation y{z) = dEg{z)/dE: 

,Jdbo{E')\ 



y{z) = 1+ r -^y{z') + ^ r dz\l + z'flMz') {- 
Jo 1 + z' Ho Jo \ 



dE' 



(B6) 



E'=(l+z')E,(z') 



Correspoinding differential equation is 
1 dy{z) 1 



y{z) dz 1 + z Hq 
The solution of Eq. ( [B7D is 



±n+^)V2M&om 



dE' 



(B7) 



E'=il+z)Egiz) 



_ dEgjz) 
^^^' ^ dE ^^^^^' ^^P 



ki>'''^'^"i 



dbojE') 
dE' 



E'=(l+z')E,{z')^ 



(B8) 



where Eg{z) is an energy at epoch z. Eq. ( |B8| ) coincides with Eq.(36) from Ref. [^ 
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